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ABSTRACT 

We employ a time-dependent synchrotron-self-Compton code for modeling contemporane- 
ous multiwavelength data of the blazar 1 ES 1218-1-30.4 The input parameters of the model 
are used to infer physical parameters of the emitting region. An acceptable fit to the data is 
obtained by taking into account a stellar emission component in the optical regime due to the 
host galaxy. The physical parameters inferred from the fit are in line with particle acceleration 
due to the Fermi mechanism providing s = 2. 1 spectra. From the properties of the host galaxy 
in the optical, we estimate the central black hole mass and thus confirm that the jet power 
injected into the emission region is in the sub-Eddington regime, as expected for BL Lacertae 
objects. 
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1 INTRODUCTION 

Among the class of active galactic nuclei (AGN), blazars are 
special in showing a spectral energy distribution (SED) that is 
strongly dominated by nonthermal emission across a wide range 
of wavelengths, from radio waves to gamma rays, and rapid, large- 
amplitude variability. Presumably, these characteristics are due to 
a relativistic jet emitted at a small angle to the line-of-sight, emit- 
ting Doppler-boosted synchrotron and inverse-Compton radiation 
and thus washing out to a variable extent the emission from the 
accretion flow and host galaxy. 

The synchrotron and inverse-Compton emission could either 
result from primary accelerated electrons, from accelerated pro- 
tons, or from secondary electrons arising in electromagnetic cas- 
cades initiated by pion and pair production (e.g. |Mannheim|1993^ . 

The high peaked BL Lacs (HBLs) show a peak in their SED in 
the X-ray regime, suggesting that an inverse-Compton peak should 
occur at correspondingly high gamma-ray energies. In fact, a large 
fraction of the known nearby HBLs have already been discovered 
with Cherenkov telescopes, such as H.E.S.S., MAGIC, and VERI- 
TAS. 

In those cases where the blazars have been detected at gamma- 
ray energies, the SED shows two bumps, one at infrared-to-X-ray 
energies, and the other at gamma ray energies. So far, the reasons 
for the variation of the peak energies are unknown, and their rela- 
tion with fundamental parameters of the central engine, the black 
hole mass, spin, and accretion rate, are far from settled. 

Models of the radiation processes behind the blazar emission 
are a key issue to improve our understanding of blazars. The diag- 
nosis of the observed spectra using radiation models allows to infer 
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the physical conditions prevailing in them and to discern their re- 
lation with fundamental parameters. However, without a prescrip- 
tion of the relativistic particle spectra, the models effectively map 
radiation onto particle spectra, and the error in the inferred phys- 
ical conditions is correspondingly large. On the other hand, using 
theoretically favoured particle spectra most strongly constrains the 
physical conditions prevailing in the sources, and this may eventu- 
ally lead to a consistent physical explanation. 

The HBL 1 ES 1218+30.4 has been discovered as a candidate 
BL Lac object on the basis of its X-ray emission and has been iden- 
tified with the X-ray source 2A 1219+30.5 ( [Wilson et j:][T979i 
[Ledden et alTl[7987] >. For the first time, 1 ES 1218+30.4 has been 
observed at VHE energies using the MAGIC telescope in January 
2005 ( [Albert et al.|200 6l. Simultaneous optical data has been pro- 
vided by the KVA telescope on La Palma. Latest TeV data have 
been provided by the VERITAS telescope ( [VERITAS CoUabora-[ 
[tion: V. A. Acciari et al."|2009| ). 

Here, we present the kinetic equation and numerical code de- 
scribing the synchrotron-self-Compton emission (Sect.|2](. The par- 
ticular emphasis lies on an accurate treatment of the Klein-Nishina- 
turnover which is important at very high gamma ray energies. Pre- 
vious works (c.f. |Bottcher & Chiang ( 2002 1 and references in there) 
provided good results in the scope of HBL modelling. Our ap- 
proach is quite similar to this models except for details in the treat- 
ment of the particular processes. External Compton effects are left 
out in order to use a minimum number of parameters, which are 
sufficient for dealing with HBLs. In Sect. |3] we apply our code to 
1 ES 1218+30.4 referring in particular to the MAGIC and VERI- 
TAS data and give a set of physical parameters for the most ac- 
ceptable fit. Finally, we discuss our results in the light of particle 
acceleration theory. 
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2 THE MODEL 



2.1.3 Compton scattering 



To model the observational data we use the well-established syn- 
chrotron self-compton (SSC) model (e.g. |Maraschi et al.|1992[ l. We 
assume a spherical, homogeneous emission region - coined blob 
- containing isotropically distributed non-thermal electrons and a 
randomly oriented magnetic field. 

Due to the presence of this magnetic field the electrons emit syn- 
chrotron radiation. The photons are then scattered oif the same 
electron population via the inverse Compton process. The resulting 
spectrum shows the typical two bump structure commonly found in 
blazars. 

In the following section the governing equations of the SSC-model 
are explained. 

2.1 Photon distribution 

To determine the time-dependent spectral energy distribution of 
blazars we solve the differential equation for the differential pho- 
ton number density, obtained from the radiative transfer equation, 
including the corresponding terms with respect to SSC model. 
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2.1.1 Synchrotron radiation 



In the following context the well known ^-approximation ( [Felten &] 
|Morrison|1966[ |Schlickeiser|2002} is applied to describe the syn- 
chrotron radiation in a convenient way. Thus the synchrotron pho- 
ton production rate Rg is given by 
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with the pitch angle averaged total power Ps emitted by a single 
electron having Lorentz factor y (^Rybicki & Lightman|1979||Blu-| 
[menthal & Gould] 1970[|Ginzburg & Syrovatskii| 19691 1 
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obtained from the pitch angle averaged critical synchrotron fre- 
quency. 

2. 1.2 Synchrotron-self-absorption 

In optically thick regimes the emitted synchrotron radiation is ab- 
sorbed by the emitting electrons itself. This is described by the syn- 
chrotron self absorption coefficient, 
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which leads to the absorption rate 
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The second main feature of the SSC model is Compton scattering 
of the synchrotron photons by the emitting electrons themselves. 
Here the full Klein-Nishina cross section is used to calculate the 
photon production rate. 
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The formula was taken from |Pe'er & Waxman| ( (2005} with minor 
corrections according to |Coppi & Blandford jl990[ l. The photon 
energies are rewritten in terms of the electrons rest mass, so that 
hv = emc- for the scattered photons and hv = e^mc^ for the target 
photons. To make use of the full Klein-Nishina cross section we 
applied the approximate inverse Compton spectrum jJones|1968| l 
of a single electron scattered off by a unit density photon field. 
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where ^" = e/(46i7^(l-e/7)) and 1/(47^) < q" < 1 . This equation 
is valid 

for 6i < e < 4ei 7-/(1 + 46i7). The corresponding ordinary Comp- 
ton spectrum is approximately given by 
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with q' = 4y^e/ei and target photon energies in the range ei/47- < 
e < 61. 

2.1.4 Photon escape 

The last term describing the evolution of the photon number density 
represents the photons escape rate. Here the photon escape time 
Tesc is given by the light crossing time 

« -, (10) 
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where Rt, is the radius of the emitting blob. The escape time is cho- 
sen to be the light crossing time of the photons. 

2.2 Electron distribution 

The time evolution of the electron distribution is described by the 
kinetic equation 
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The synchrotron loss is given by ys = Ps(y)/inc- with the syn- 
chrotron power Ps(y) (cf. l[3|). T^^ce- = t]Rblc describes the elec- 
trons escaping from the emission region, where rj is an empirical 
factor. The inverse compton losses yic including the full Klein- 
Nishina cross-section are adopted following |SchlickeiserH200^ . 
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Tmax B/G K/cm s R/cm S rj 

5.0x10^ 2.1 0.04 0.4x 10"' 3 x lO'^ 80 10 

Table 1. Best fit parameters for the SSC modelling the SED of 1 ES 
1218+30.4 



with 
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As an injection function Qi„j we use a power law combined 
with an exponential cut-off: 



Qmjiy) = Ky exp 
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with Lorentz factor y, the normalisation constant K, cut-off energy 
ymax and the spectral index s. Using a injection function constant in 
time one yields a equilibrium solution for n^- (y) being constant in 
time. 



2.3 Numerics 

To obtain a model SED using the SSC formalism we solve the cou- 
pled equations eq. l|TJ and eq. numerically in our code frame- 
work. With respect to stability issues we use the Crank-Nicholson 
scheme ( |Press|2002| l to compute the synchrotron part of the right 
hand side of both electrons equation and photons equation. The 
code was tested carefully and stands the challenge of computing 
the equations in a range of 20 orders of magnitude. All single 
effects (synchrotron radiation/losses, Compton scattering/losses) 
have been cross checked with analytical solution and approxima- 
tion as well with numerical integration cross checks using Math- 
ematica. Also comparison with existing codes were done success- 
fully (as long as the models themselves were comparable). 



3 RESULTS 

In Fig[T|we show the results of the application of the code to the 
data of the HBL 1 ES 1218-1-30.4. Here we present a SSC model 
curve which fits the data in the X-ray and VHE regimes. The corre- 
sponding parameters are listed in Tab.[T] With this values we end up 
in an equilibrium state with a cooling break energy of the electrons 
distribution y break = 6 x 10*. 

The VHE data shown here have been discovered by the 
MAGIC telescope ( {Albert et al.|2006^ . Lately the VERITAS tele- 
scope could confirm this detection. As shown in |VERITAS Collab^ 
loration: V. A. Acciari et al. | P009[ | the measured flux matches each 
other in the overlapping energy regime. 

The data in the X-ray regime have been taken by SWIFT between 
March and December 2005 ( [Tramacere et al.||2007[ l. Another set 
of X-ray data were obtained by the BeppoSAX experiment in 1999 
( [Donato et al.|2005] l. It is remarkable, that in spite of a lag of 6 years 
between these observations there is no difference in the X-ray flux. 
Together with the constant TeV flux level this is a strong argument 
for an almost constant background of non-thermal electrons in a 
constant magnetic field. 



The KVA data point shown in Fig. [T] was obtained simulta- 
neously to the MAGIC observation of 1 ES 1218-1-30.4. One can 
also see that the SSC model is not able to fit this data point. This 
discrepancy can be resolved by taking the NED and 2MASS data 
( |Chen et al.|2005| l surrounding the optical KVA point into account. 
This set of data points has been modeled by a simple blackbody 
spectral distribution with a temperature T = 4500K and a ra- 
dius Ru, = 2.85 X 10'*cm. Considering average sun-like stars we 
yield an estimate of the number of stars responsible for this black- 
body radiation given by R^jJE^ = 1.7 x 10". This number to- 
gether with the temperature T gives a hint for the host galaxy of 
1 ES 1218H-30.4 being the origin of this feature in the SED not be- 
ing modelled by the SSC approach. With this blackbody approx- 
imation we get an estimate for the central black hole mass of the 
AGN using the correlation given by |Kormendy & GebhardtlpOOl) 
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Here we get Mg^ = 5.6 x 10** and a corresponding Eddington lu- 
minosity Lejj = 7.3 X lO"**. Lacking observations of the radial de- 
pendence of the surface brightness of the host galaxy, the accuracy 
of these estimates is not better than a factor of a few. Adopting a 
bulk Lorentz factor F = 57 , which is consistent with <5 = 80 for 
reasonable angles of the jet axis to the line-of-sight, the injected lu- 
minosity in the AGN frame is L,„^ = Aj'inR^T- j dyymc^ Qiy) = 
7.3 X 10"*^^ erg/s. The resulting very low Eddington ratio of the order 
of 0.001 is in line with the results of population studies of BL Lac 
objects, for which [T7eves et al^ ( |2002| l find 0.01. 



4 DISCUSSION 

In this paper, we presented SSC model fits to the contemporaneous 
data of 1 ES 1218+30.4 in the X-ray and VHE range. In the optical 
regime a simple blackbody spectrum has been applied to complete 
the model SED. In a similar way, |Katarzyilski et al.| ( |200T| l modeled 
the optical data of Markarian 501 by applying a standard model for 
the elliptical host galaxy ( [Nils son et a l. 1999 ). Our simple approach 
to the background radiation of the host galaxy actually suffices, the 
more sophisticated approach of Katarzyriski et al._ ( 2001 1 might be 
of use if more data for the host galaxy would be available, but at this 
high redshift details are not accessible. [Celotti & GhiselliniH2008] > 
have recently studied 1 ES 1218+30.4 , assuming it is in a flaring 
state. 

Observational results for 1 ES 1218+30.4 have also been dis- 
cussed in [Sato et aT] ^2008^ with special regard to the variability 
of the source. The authors concluded that the source must have a 
very hard electron distribution with power law slope s = 1.7. We 
disagree with this result, as we were able to show that there ex- 
ist model parameters which are well in line with relativistic shock 
acceleration theory ( [Ellison et al.||199^ , although harder spectra 
can be imagined for more extreme sources ( [Vainio & Schlickeiser| 
[T998l l. 

The major difference of our spectrum compared to |Sato et al.| 
l |2008| l is that we assumed the optical regime to be dominated by 
the host galaxy, approximately described by a blackbody spectrum. 
Therefore, the need for extreme electron spectra could be relaxed. 
Considering that the SWIFT, MAGIC and VERITAS data used here 
do not show strong flaring features, we have modeled the SED as 
steady-state SSC emission with our time-dependent code, obtaining 
physical parameters of the emission region. These parameters lie 
well in the range found with SSC models for other HBLs. The small 
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Figure 1. Overall SED of 1 ES 1218+30.4 . Green crosses show the MAGIC data ^Albert et al.|2006) , grey triangles the VERITAS data I VE RITAS Collab-| 
[oration: V. A. Acciari et al. [2009^ . The TeV data have been de-absorbed by applying the EBL correction model of [Kneiske et al. (20041. The filled cyan 
box represents the simultaneously obtained KVA data point. In the X-ray range SWIFT (Tramacere et al. 2007i and BeppoSAX (Donato et al.]2005^ data are 
plotted as blue stars and red crosses. The 2MASS data (red open triangles, (Chen et al. 2005 together with the NED data represent the peak in the optical 
range fitted by a thermal spectrum (dashed green line). The red solid line depicts the SSC model curve. 



magnetic field value differs slightly from the normally used values 
of about 0.1 - 0.2 G in SSC models, but in contrast to competing 
hadronic models this value is still reasonable. 

Additionally, the central black hole mass could be estimated 
from the host galaxy properties, demonstrating that the emission 
region is consistent with a sub-Eddington jet as generally expected 
for BL Lac-type sources.The growing, but still marginal, discrep- 
ancy of our model SED and the VHE spectra at highest energies, 
if taking the expected gamma ray attenuation due to pair produc- 
tion into account, is a concern. It is also found in other studies of 
TeV blazars. This trend could indicate an insufficiency of the SSC 
model approach, a weaker than expected gamma ray attenuation, 
or an incomplete understanding of the energy determination of air 
showers from their Cherenkov emission. 
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